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Abstract 



The perturbative expression of the running strong couphng constant as((5^) 
has an unphysical singularity for = Aq(jj-,. Various modification have been 
proposed for the infrared region. The effect of some of such proposals on 
the quark-antiquark spectrum is tested on a Bethe-Salpeter (second order) 
formalism which was successfully applied in previous papers to an overall 
evaluation of the spectrum in the light-light, light-heavy and heavy-heavy 
sectors (the only serious discrepancy with data being for the light pseudoscalar 
meson masses). In this paper only the cc, bb and qq (q = u or d) cases are 
considered and fine structure is neglected. It is found that in the bb and cc 
cases the results are little sensitive to the specific choice. In the light-light 
case the Dokshitzer et al. prescription is again essentially equivalent to the 
truncation prescription used in the previous calculation and it is consistent 
with the same a priori fixing of the quark light masses on the typical current 
values mu = = 10 MeV (only the pion mass resulting completely out 
of scale of about 500 MeV). With the Shirkov-Solovtsov prescription, on the 
contrary, a reasonable agreement with the data is obtained only at the price of 
using a phenomenological momentum dependent effective mass for the quark. 
The use of such an effective mass should amount to a correction of the free 
quark propagator. It is remarkable that this has also the effect of bringing 
the pion mass in the correct range. 



PACS: 12.38.Aw, 11. 10. St, 12.38.Lg, 12.39.Ki 
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I. INTRODUCTION 



In perturbation theory the running couphng constant in QCD is usually written up to 
one loop as 

4:71 



or also up to two loops 



/5oln(QVA2 



1 + 



2A ln(ln(gVA^)) 
/?o' ln(gVA2) 



(2) 



Q being the relevant energy scale, Po = 11 — |nf, /?i = 51 — yrif and rif the number of 
flavors with masses smaller than Q. 

Such expressions have been largely tested in the large Q processes and are normally used 
to relate data obtained at different Q using the appropriate number of "active" flavors nf 
and different values of A in the ranges between the various quark thresholds. 

Both expressions become singular and completely inadequate as approaches A^. 
Therefore they must be somewhat modified in the infrared region. 

Various proposals have been done in this direction. The most naive assumption consists 
in cutting the curve (|I|) at a certain maximum value as(0) = as to be treated as a mere 
phenomenological parameter (truncation prescription). Alternatively, on the basis of general 
analyticity arguments, Shirkov and Solovtsov [|l| replace (|l]) with 

This remains regular for = h? and has a finite A independent limit as(0) = 47r//3o, for 
0. Finally, inspired also by phenomenological reasons, Dokshitzer et al. [Q write 

= a's{Q% (4) 

where a^lQ"^) is the perturbative running coupling constant as given by Eq. (Q) and V = 
d / d{\n.{Q'^ / Is?')) is a derivative acting on al^Q"^). The various curves are reported in Fig. |l|. 

The above modified expressions have been applied to study various effects in which 
infrared behavior turns out to be important. Electron-positron annihilation into hadrons, 
r-lepton decay, lepton-hadron deep inelastic scattering, jet shapes, pion form-factors etc. 
are of this type. 

In the quark-antiquark bound state problem the variable can be identified with the 
squared momentum transfer = (k — k')^ and formally the use of a running coupling 
constant amounts to include higher order terms in the perturbative part of the potential 
or the Bethe-Salpeter kernel. In this case all values of are involved and an infrared 
regularization becomes essential. Furthermore (Q^) ranges typically between (1 GeV)^ and 
(0.1 GeV)^ for different quark masses and internal excitations and values of smaller than 
A^ can be important. The specific infrared behavior is therefore expected to affect the 
spectrum and other properties of mesons. 

The purpose of this paper is to test such kind of effects in a particular formalism we have 
developed and used in previous papers. 
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II. FORMALISM 



In reference [Q] we have obtained a good reproduction of the entire meson spectrum in 
terms of only four adjustable parameters, by solving numerically the eigenvalue equation for 
the squared mass operator 



= + U 



(5) 



or the mass operator 



M = Mn + V, 



(6) 



where Mq = wi + 1102 = \J rrii + + y m2 + is the kinetic term and U and V are com- 
plicated momentum dependent potentials. Up to the first order in the running coupling 
constant as(Q^) and in terms of the string tension a, the "quadratic potential" U is given 
by 



(k|f/|k') 
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with 



J"^''(r,q, gio,g2o) 
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k-k'. 



and a'; = 7°7^ = \ e^^^^lj, 7]"], Qjo = {w, + w'^)/2, j = 1, 2. Q 
q = (k + k') /2, q\ = {6''^ - f^r^)q^^ v = r/r. 

The above expression was obtained by reducting of a Bethe-Salpeter like equation which 
was obtained in reference from first principle QCD under the only assumption that the 
logarithm of the Wilson loop correlator W could be written as the sum of its perturbative 
expression and an area term 



ilnW = i(lniy)pert + aS 



(9) 



(advantage was taken in the derivation of an appropriate Feynman-Schwinger representation 
of the quark propagator in an external field). 
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An expression for (k|l^|k') can be obtained by a direct comparison of Eq. with 
Eq. (P). Neglecting terms in V"^ (what is consistent with the other approximations), this 
amounts simply to change the kinematical factor in front of the right-hand side of Eq. (^. 
Properly one should divide (k|[/|k') by Wi+W2+w[+W2, practically the simpler replacement 
{wi+w2){wj^+w2) _^ — — 1 — essentially equivalent. 

The interest of the more conventional Eq. (|^) is that it makes more immediate a compar- 
ison with ordinary potential approaches and the consideration of the non-relativistic limit. 
In particular V coincides with the Cornell potential (k|V|k') = (k|(— |^ + crr)|k') in the 
static limit and with the potential obtained in P] when the first relativistic corrections are 
included. In this paper, however, we shall refer only to Eq. (^, as more directly related to 
the original B-S equation. 

The method used in consists in solving first the eigenvalue equation for M in the 
static limit of V by the Rayleigh-Ritz method (using an harmonic oscillator basis); then in 
evaluating the quantity (M^) (or (M)) for the resulting zero order eigenfunctions. 

Actually in we have neglected the complicated fine structure spin dependent terms 
occurring in (|^ and taken into account only the hyperfine term in - cr^- We have used 
the expression (|l]) for the running coupling constant with Uf = 4 and A = 0.2 GeV frozen 
at as = 0.35; we have also taken a = 0.2 GeV^, mu = md = 10 MeV Q, rric = 1.394 GeV, 
mb = 4.763 GeV. The quantities A and = were fixed a priori from high energy data; 
Ug, cr, rric and mb were adjusted on the ground cc and bb states, the cc hyperfine splitting 
and the Regge trajectory slope. 

In this paper, to test the sensitivity of the results to the infrared behavior of as(Q^), we 
have performed the same calculation in the bb, cc and qq {q = u or d) case using Eqs. (H) 
and (^), with an appropriate redefinition of the adjustable parameters. 

III. RESULTS 

In tables |I|, 0, and|ITT|we give the bb, cc, and qq quarkonium masses respectively obtained 
for the different running coupling constant prescriptions. In column (a) we report the results 
obtained in ref. for the truncated as(Q^), in column (b) those obtained by means of Eq. 
(0) proposed by Shirkov-Solovtsov and in column (c) those obtained by means of the q;s(Q^) 
of Eq. (^) proposed by Dokshitzer et al. 

In columns (b) and (c) we used the same values nf = 4, A = 0.2 GeV and m-a = 
rrid = 10 MeV as in 0], but we have slightly redefined the adjustable parameters taking 
cr = 0.18 GeV^ in both cases and = 1.545 GeV and mb = 4.898 GeV for prescription (^ 
(column (b)), rric = 1.383 GeV and mb = 4.7605 GeV for prescription @ (column (c)). 

Notice that, in spite of the reduced number of adjustable parameters, the spectra of 
bottonium and charmonium are not essentially modified by the new choice for as(Q^), with 
perhaps the exceptions of the highest cc states that are lower in the Dokshitzer et al. case. 
This indicate little sensitivity of such spectra to the infrared behavior of q;s(Q^)- 



Notice that the results are very little sensitive to the precise values of and md if these are 
small. 
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The situation is completely different for the light-light spectrum of table |TT|. While in 
front of the experimental and theoretical uncertainties columns (a) and (c) can be considered 
not really distinguishable, the values reported in column (b) are definitely systematically 
too low (in particular (M^) < for vr-meson). 

Notice however that in the above reported calculations we have used for mu and rrid the 
current mass value of 10 MeV. This amounts to assume the difference between the current 
and the constituent masses to be essentially related to kinematical relativistic correction (cf. 
0) or, what is the same, that the free quark propagator is a good approximation for the 
complete one in the B-S equation. The inability of the formalism to reproduce a reasonable 
value for the vr mass and all experience gained by the chiral symmetry problematic (see in 
particular ^ and reference herein) suggests that this should not be the case for the light-light 
systems. 

For this reason we have repeated the calculation for choice (^ with various constituent 
values for the light quark masses. Two sets of results are reported in columns (d) and (e) 
of table 1^. Notice that for = = 0.30 GeV the situation is again very similar to 
those of column (a), table |IT|. Notice also however that as mu.d increases the bound state 
masses uniformly increase and that for low value of mu,d the lowest bound state masses can 
be made to agree fairly well with the data, for high value the same occurs for higher states. 
This could suggest the use of a kind of running constituent mass. 

In column (h) the results are reported for the squared effective mass 

mis = 0.11 k- 0.025 + 0.265 (10) 

k denoting the quark momentum in the center of mass frame. In Eq. (|IU]) the coefficients 
are chosen in order to obtain m^s = 0.22, 0.28, 0.35 GeV for k^ = 0.26, 0.41, 0.58 GeV^ 
approximately corresponding to the (A;^) values for the IS, 2D and IG states respectively. 
As it can be seen the agreement with the data is much improved in this way and finally even 
a reasonable value for the vr mass is obtained. 

Notice that the use of an effective running mass is in agreement with the general perspec- 
tives of combining Dyson-Schwinger equation with B-S equation. Obviously a fine tuning 
of the coefficients in Eq. (|I^) to further improve the results would be meaningless in the 
present context, due even to the approximation used {e.g. the triplet-singlet splitting has 
been essentially evaluated perturbalively) . 

IV. CONCLUSIONS AND DISCUSSION 

In conclusion the heavy quarkonium spectrum does not seem to be very sensitive to the 
specific infrared behavior of the running constant as it could be expected on general ground. 

For the light quarkonium Dokshitzer et al. prescription (^) does not essentially change 
the results in comparison with the truncation assumption adopted in |^, in spite of the very 
different appearance of the corresponding curve in Fig. [l| 

On the contrary the situation changes drastically for prescription (|^). 

Actually such prescription would be definitely ruled out if we insisted in using a current 
mass for the light quarks. If however we give to m^^^ a constituent value the results are 
again similar to those obtained with the truncated as(Q^) and can be strongly improved if 
we use the running effective mass (0). 
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To understand better the meaning of Eq. (0), notice that, in the context of the second 
order formahsm developed in ref . ||^ , the free quark propag ator occurring in the BS 

equation is i/ip^ — rn?), where we have to set p = {Mb/2 ± ko, ±k) (the upper and lower 
signs referring to the quark and the antiquark respectively) in the CM. frame . However if, 
consistently with the other approximations made, we neglect the spin dependent terms, the 
full quark propagator i?2(p) can be written as i/{p'^ — + T{p)). Then, recalling that the 
instantaneous approximation consists in setting fco = in the BS kernel, in the same order 
of ideas we can replace the slowly varying expression T{po, |p|) by T{Mb/2, |k|). Eventually 
we obtain the operator as given by Eq. (^) but with mi = m2 replaced by 

m,\(|k|)=m2-r(MB/2,|k|). (11) 



Eq. (0) corresponds to a parametrization of right hand side of Eq. (11). Obviously 
in principle T{p) should be obtained by actually solving the DS equation for H2{p) (ref. 
[3b,4b,4c]), but this is a complicated task that we reserve to a forthcoming paper. 
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FIGURES 




FIG. 1. Running coupling constant as{Q) on logarithmic scale. Truncation prescription (full 
line), Shirkov-Solovtsov prescription (dashed line), Dokshitzer et al. prescription (dot-dashed line). 
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TABLES 



TABLE 1. bb. Uf = A, A = 0.2 GeV. (a) mt = 4.763 GeV, a = 0.2 GeV^ a^iO) = 0.35, 
Ref. [3]. (b) mb = 4.898 GeV, a = 0.18 GeV^ Shirkov-Solovtsov a^iQ'^). (c) mb = 4.7605 GeV, 
a = 0.18 GeV^, Dokshitzer et al. a^iQ^). 



States 




exp. 




fa) 










(GeV) 




fGeV) 


fGeVl 


(GeV) 


1 bo 








y.o ( 4 


y.o / 4 


O Q7F^ 


1 3q 
1 bi 


I(li) 


9.46037 ± 0.00021 


y.4DU 


n /I i^n 
y.4DU 


y.4DU 


9 Ici^ 
Z Oq 










Q QS8 
y.yoo 


y.yoo 


3c! 




10.02330 ± 0.00031 


lU.UiU 


1 n n9Q 
iU.Uzo 


1 n ni v 
iU.Ui ( 










lU.OZZ 


1 n ^A9 


1 n ^98 

lU.OZO 


3q, 
o oi 


1 (^d^j 


10.3553 ± 0.0005 


iU. o4o 


iU.oDo 


iU.oOZ 


4 lc;„ 

4t Do 








1 n 


1 n fil 8 




/I 3q 
4 oi 


1(4^) 


10.5800 ± 0.0035 


lU.DzU 


1 n ftQQ 
iu.Doy 


1 n fji 9 
iU.DiZ 


DO 








iU.OO / 


1U.C504 




c 3q 


1 (^iUoou ) 


10.865 ± 0.008 


iU.OO ( 


iU.o / / 


iU.oo4 


e^So 








11.060 


11.070 


11.026 




T(11020) 


11.019 ±0.008 


11.079 


11.089 


11.044 


i^Pi 








9.908 


9.918 


9.914 


l^P2 


Xb2(lP) 


9.9132 ± 0.0006 ' 










l3Pi 




9.8919 ± 0.0007 


> 9.900 


9.908 


9.920 


9.917 


l3Po 


X6o(lP) 


9.8598 ± 0.0013 










2^Pi 








10.260 


10.279 


10.269 


2 3P2 


x;.2(2P) 


10.2685 ± 0.0004 ^ 










2 3Pi 


X6i(2P) 


10.2552 ± 0.0005 


> 10.260 


10.260 


10.280 


10.271 




\w-)(2-P) 


10.2321 ± 0.0006 _ 
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TABLE 11. 


cc. n/ = 4, A 


= 0.2 GeV. (a) 


mc = 1.394 GeV, a 


— 2 GoV^ 


n (0) — Q 35 


Ref. [3]. (b) mc = 1.545 GeV, 


cr = 0.18 GeV^ Shirkov-Solovtsov a. 


■ (0'^) (c) 


= 1 383 GeV 


a = 18 GeV^ 


Dokshitzer pf: o.l 


. «s(Q'). 








States 




exp. 


(a) 


(b) 


(c) 






(MeV) 


(MeV) 


(MeV) 


(MeV) 


l^So 




2979.8 ± 2.1 


2982 


2977 


2982 


l^Si 


J/^(15) 


3096.88 ± 0.04 3097 


3097 


3097 


lASS 




117 


115 


119 


116 


2^So 


r/c(25) 


3594 ± 5 


3575 


3606 


3573 




^^(25) 


3686.00 ± 0.09 3642 


3670 


3636 


2 ASS 




92 


67 


64 


63 


S^So 






3974 


4005 


3950 


S^Si 


V'(4040) 


4040 ± 10 


4025 


4054 


3998 


4IS0 






4298 


4323 


4252 


43Si 


^^(4415) 


4415 ± 6 


4341 


4364 


4291 


l^Pi 






3529 


3556 


3528 


13P2 


Xc2(lP) 


3556.17 ±0.13 1 








l3Pi 


Xcl(lP) 


3510.53 ±0.12 


} 3525 3530 


3561 


3531 


I'Po 


Xco(lP) 


3415.1 ± 1.0 J 








2iPi 






3925 


3954 


3904 


2 3p 






3927 


3958 


3906 


1^D2 






3813 


3853 


3811 


1 3f»„ 
i JJ3 












l^Da 


V'(3836) 


3836 ±13 


1 3813 


3854 


3811 


l^Di 


^^(3770) 


3769.9 ± 2.5 J 








2ID2 






4149 


4183 


4121 


2^03 












2^02 






1 4149 


4184 


4121 


2^01 


V'(4160) 
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TABLE III. qq. m^.d = 0.01 GeV, n/ = 4. (a) as(0) = 0.35, Ref. [3], a = 0.2 GeV^, A = 0.2 
GeV. (b) Shirkov-Solovtsov as(Q^), = 0.18 GeV^ A = 0.2 GeV. (c) Dokshitzer et al. a^iQ'^), 
a = 0.18 GeV2, A = 0.2 GeV. 



States 




exp. 
(MeV) 






(a) 
(MeV) 


(b) 
(MeV) 


(c) 
(MeV) 


l^Sn 


I"! 


134.9764 ± 0.0006 1 


479 




575 






139.56995 ± 0.00035 J 




l^Si 


p(770) 


768.5 ± 0.6 




846 


423 


904 


lASS 




630 






367 




329 


2^So 


7r(1300) 


1300 ± 100 


1326 


952 


1338 


23Si 


p(1450) 


1465 zb 25 




1461 


1128 


1459 


2 ASS 




165 






135 


176 


121 


S^So 


7r(1800) 


1795 ± 10 




1815 


1485 


1793 


33Si 


p(2150) 


2149 ± 17 




1916 


1600 


1889 


3 ASS 




354 






101 


115 


96 


l^Pi 


6i(1235) 


1231 ± 10 










1=^P2 
l=^Pl 


02(1320) 
ai(1260) 


1318.1 ±0.7 ' 
1230 ± 40 


> 1303 


1333 


1045 


1365 


l^'Po 


ao(1450) 


1450 ± 40 












l^Ds 


7r2(1670) 


1670 ± 20 










l^Ds 

1^D2 


P3(1690) 


1691.1 ±5 




> 


1701 


1444 


1715 


l^Di 


p(1700) 


1700 ± 20 












I'Fg 
















1=^F4 
l^Fg 


04(2040) 
X(2000) 


2037 ± 26 


> 




1990 


1743 


1985 


1=^F2 






> 










1^G4 
















1=^G5 
13G4 


P5(2350) 


2330 ± 35 


> 




2238 


1994 


2214 


l^Gs 


P3(2250) 




> 










l^Hs 
















l^He 
l^Hs 


06(2450) 


2450 ± 130 


> 


2460 


2215 


2416 


13H4 
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TABLE IV. qq. Uf = 4. Shirkov-Solovtsov as(Q^). a = 0.18 GeV^. A = 0.2 GeV. (d) 
mu,d = 0.22 GeV. (e) mu,d = 0.30 GeV. (f) ^ = 0.11 - 0.025 fc^ + 0.265 k'^. 



States 


(McV) 


exp. 






(d) 


(e) 


(f) 


I'So 




134.9764 ± 0.0006 1 
139.56995 ± 0.00035 J 


26 


473 


124 


1 ^Si 
J. j\ 




768.5 ± 0.6 




725 


868 


737 


A SS 




630 






699 


394 


613 


2^80 


7r(1300) 


1300 ± 100 




1190 


1326 


1401 






1465 ± 25 




1 ^44 


1 468 




A SS 




165 






154 


142 


107 


S^Sn 


7rfl800) 


1795 ± 10 




1688 


1806 


1993 


S^Si 


p(2150) 


2149 ± 17 




1788 


1900 


2063 


A SS 




354 






1 0(1 


Q 1 


7(1 


J- 1 


6i fl235) 


1231 ± 10 










l^Po 
l3Pi 


09 f 1320) 
ai(1260) 


1318.1 ±0.7 ' 
1230 ± 40 


> 1303 


1243 


1364 


1319 


l^'Po 


ao(U50) 


1450 ± 40 












1 ^Do 


TTofieTO") 


1670 ± 20 










I^Dq 

1^D2 




1691.1 ±5 




> 


1603 


1715 


1741 


1 ^Di 


n(^ 7nn'i 


1700 ± 20 












1 ^F-i 
















l^F/i 

J. J. 4 
l^Fg 

1 ^Fo 


aA (2040^) 
X(2000) 


2037 ± 26 


> 


> 


1881 


1979 


2043 


1^G4 
















l^Gs 

13G4 


/95(2350) 


2330 ± 35 


> 


> 


2118 


2209 


2319 


l^'Gg 


P3(2250) 














l^Hs 
















I'He 
l^Hs 


06(2450) 


2450 ± 130 


> 


2329 


2415 


2569 


13H4 

















11 



REFERENCES 



[1] D.V. Shirkov and I.L. Solovtsov, Phys. Rev. Lett. 79 (1997) 1209; 
Theor. Math. Phys. 120 (1999) 1220; 

N.G.Stefanis, W. Schroers, H.-Ch. Kim, Phys. Lett. B449 (1999) 299. 
[2] Yu.L. Dokshitzer, A. Lucenti, G. Marchesini, G.P. Salam, JHEP 9805 (1998) 003; 

Yu.L. Dokshitzer, in 29*^^ International conference on High-Energy Physics (ICHE 98), 
Vancouver, Canada, A. Astbury, D. Axen, J Robinson Eds. World Scientific (1999), [tiep- 
ph/981225^ ; 



Yu.L. Dokshitzer, V.A. Khoze, S.I. Troyan, Phys. Rev. D53 (1996) 89. 
[3] M. Baldicchi and G.M. Prosperi, Phys. Rev. D62 (2000) 114024; 
Fizika B 8 (1999) 2, 251; 

M. Baldicchi, in '^QCD: Perturbative or Nonperturbative ?" Pag. 325, L. S. Ferreira, P. 
Nogueira, J. I.Silva-Marcos Eds. World Scientific (2001), |hep-ph/9911268 . 



[4] N. Brambilla, E. Montaldi and G.M. Prosperi, Phys. Rev. D54 (1996) 3506; 
G.M. Prosperi, (1997), |E^h/970904q ; 

G.M. Prosperi, in Problems of Quantum Theory of Fields, Pag. 381, B.M. Barbashov, 
G.V. Efimov, A.V. Efremov Eds. JINR Dubna 1999, |hep-ph/9906237| . 
[5] A. Barchielh, E. Montaldi and G.M. Prosperi, Nucl. Phys. B296 (1988) 625; 
Erratum-ibid. B303 (1988) 752; 

A. Barchielh, N. Brambilla and G.M. Prosperi, II Nuovo Cimento 103A (1990) 59; 
N. BrambiUa, P. Consoli and G.M. Prosperi, Phys. Rev. D50 (1994) 5878. 
[6] M.B. Hecht, C.D. Roberts, S.M. Schmidt, Phys. Rev. C63 (2001) 025213, 
ANL-PHY-9742-TH-2000, |nucl-th/0011005t 



ANL-PHY-9706-TH-2000, |nucl-th/0010"024 



C.D. Roberts, ANL-PHY-9631-TH-2000, |nucl-th/0007UM 



12 



